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Abstract 

If the scintillator response to a hadronic shower in a semi-infinite uniform 
calorimeter structure is S relative to the electronic response, then S/E = 
[fern + (1 ~ /eTn)(^/e)], where E is the incident hadron energy, /em is the elec- 
tronic shower fraction, and h/e is the hadron/elcctron response ratio. If there is 
also a simultaneous readout with a different h/e, say a Cherenkov signal C, then 
a linear combination of the two signals provides an estimator of E that is propor- 
tional to the incident energy and whose distribution is nearly Gaussian — even 
though the S and C distributions are non-linear in E, wide, and skewed. Since 
c/3 ■ an estimator of fem is also obtained, it is no longer a stochastic variable. Much 

^> of the remaining resolution variance is due to sampling fluctuations. These can 

be avoided in a homogeneous calorimeter. The energy resolution depends upon 
the contrast in h/e between the two channels, h/e is small in the Cherenkov 
channel. Mechanisms that increase h/e in sampling calorimeters with organic 
scintillator readout are not available in a homogeneous inorganic scintillator 
calorimeter. The h/e contrast is very likely too small to provide the needed 
energy resolution. 
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1. Introduction 

A homogeneous dual-readout hadron calorimeter has been suggested for pos- 
sible use at a future linear collider The machine will probably be an e'^e~ 
collider and in some concepts will have a long bunch spacing (0(100 ns)). so 
that detectors with time constants in this range might be used. Discrimination 
between the Cherenkov signal (C) and scintillator (S) optical signals is expected 
to use a combination of timing, color, and, possibly, Cherenkov light direction 
and polarization. 

In practice corrections must be made for cracks, leakage, and light collection 
variations, and the structure usually varies with depth. For the purposes of 
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this analysis we shall assume that the corrections have been made properly, 
and consider a semi-infinite calorimeter with uniform structure that is either 
fine-sampling or homogeneous. 

In each high-energy interaction in a hadronic cascade an average of 1/4 of 
the energy is carried away by tt^'s 0. These immediately decay to 7's which 
initiate electromagnetic (EM) showers. This occurs many times, with the result 
that a large fraction fem of the incident energy joins the EM shower. The mean, 
if em), is ~ 0.5 for 100-150 GeV incident pious. It increases slowly with incident 
energy E, asymptotically approaching unity. 

The hadronic response S to an incident hadron with energy E (calibrated 
to electron response) is 

S^E[U„ + il-fern){h/e)]. (1) 

The EM energy deposit is detected with relative efficiency e, and the hadronic 
signal with relative efficiency h. Both vary from event to event. In part because 
of low multiplicities in the initial hadronic interactions, the variance of h is 
much larger than the variance of e. It makes sense to treat h/e as a stochastic 
variable. To the extent that the variance of h dominates, the distribution of 
the conventional e/h is not useful. In Sec. [5] we treat the distribution of h/e as 
Gaussian. 

Most energy deposit is by very low-energy electrons and charged hadrons. 
Because so many generations are involved in a high-energy cascade, the hadron 
spectra are essentially independent of the cascade's origin except for overall nor- 
malization. This "universal spectrum" concept is discussed in detail in Ref. 
It is because of this feature that (h/e) is a robust quantity, independent of 
energy and incident hadron species. 

The energy-independent (h/e) does depend upon calorimeter composition 
and structure, as well as the readout — for example, an organic scintillator read- 
out is sensitive to the otherwise-invisible neutron content of the cascade while 
a Cherenkov readout is relatively blind to the hadronic content, (fem) can be 
found by fitting the average 7r~/e response as a function of test-beam energy 
with an appropriate (fem) parameterization such as a power law in energy (^jlj 

Usually {h/e) is less than unity, since the EM contribution is detected with 
greater efficiency than the hadronic energy deposition. If (h/e) is not unity, then 
the broad, skewed fem probability distribution function (p.d.f.) significantly 
degrades and skews the energy resolution, resulting in the familiar wide, non- 
Gaussian energy distributions. The response is not linear with energy because 
of the energy dependence of (fem)- If fem could be measured for each event, 
then the response as given in Eq.[T] could be corrected to the actual energy with 
a nearly Gaussian distribution and a mean proportional to the energy. 



^Technically, a power-law fit finds a = (1 — (/i/e))i?Q^'". Since 1 — m is small and the scale 
energy Eg is close to 1 GeV for pion-induced cascades, the distinction is minor: {h/e) fa 1 — a. 
A similar distinction occurs when other paramcterizations are used, (h/e) itself cannot be 
isolated. 
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The importance of measuring the EM content on an event-by-event basis 
was reahzed as early as 1980, although how to use the information was not so 
clear. There was even a (stillborn) dual readout test by A. Erwin (BNL) using 
scintillator and radiator plates Q. 

EM showers result in large local energy deposit; with sufficient readout seg- 
mentation this "lumpincss" provides a measure of fem- Weighting this part 
differently than the remaining signal might improve resolution. This approach 
was used with some success by the WAl collaboration but has been less 
successful elsewhere, e.g., for the ATLAS central barrel calorimeter 

In a 1983 summer school review of high-energy calorimetry, P. Mockett 
stressed the importance of measuring the fractional EM content of the shower. 
He speculated that one could use two sampling media, an electron-sensitive de- 
tector (Cherenkov) and an ionization sensitive detector (scintillator). He also 
imagined taking advantage of the fast Cherenkov pulse and slow scintillation 
signal in a heavy inorganic scintillator. Both suggestions were prophetic Q. 

Such a separation was actually made by Thcodosiou et al. [7[ in 1984, using 
the time structure of pulses observed in scintillating glass. He thought the 
technique might permit electron/hadron separation or even help with particle 
identification. Winn later suggested using color in addition to timing to make 
the separation Q . There must have been considerable speculation about dual- 
readout calorimetry, but only Theodosiou et al. took this speculation into the 
laboratory. 

Part of the problem was that the physics of energy deposition had not yet 
been elucidated, or at least widely understood. This came in the late 1980s 
with the work of Fabjan et al. Q, Wigmans Briickmann et al. 11 1, Drews 



et al. jl2l |. and others, but a key element was the energy deposition inventories 
produced by the very detailed simulations of Gabriel and his collaborators at 
Oak Ridge as early as 1974 (ilH. 

Much of the hadronic energy resolution problem was related to the large frac- 
tion of missing energy in the hadronic sector, due to nuclear dissociation, nuclear 
recoil, residual nuclear excitation, ^ and v escape, and (unobserved) neutrons. 
Scintillator response to highly ionizing charged particles is non-linear, resulting 
in significantly more lost signal. For a time it was thought that ionization by U 
fission products could make up some of the lost energy [^, l3j lil, but non-linear 



scintillator response to the highly ionizing fragments negated most of the gains. 

In a sampling calorimeter, only a small fraction of the energy is deposited 
in the sensors (quartz or scintillator) , and fluctuations in this fraction are more 
important than intrinsic fluctuations in the hadronic signal. These dominate, 
once fem is removed. The sampling fluctuations are avoided in a homogeneous 
calorimeter. The possibly long bunch spacing at a future linear collider opens the 
door to a homogeneous dual-readout dense crystal or glass calorimeter, where a 
fast, blue, Cherenkov pulse might be separated from a slower, redder scintillation 
signal. Crystal studies are being successfully explored by Akchurin, et al. [lil - 
|22| . but with only speculative mention of dual- readout hadron calorimeters. A 
feasibility study is part of a new proposal Q ■ 

Akchurin et al. have demonstrated signal separation that would be adequate 
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for recovering energy estimators that are linear in the corrected energy and have 
a nearly Gaussian distribution. However, I am concerned that the energy res- 
olution would not be adequate. In this paper I explore the likely resolution as 
a function of energy and {h/e\s) using resolution contributions based on pub- 
lished crystal, glass, and sampling calorimeter performance. Simple, transparent 
Monte Carlo simulations (MC's) are used by choice, to make the physics more 
transparent than if a sophisticated MC such as GEANT4 were used. The p.d.f. 
of fern IS approximated with some care, while other resolution contributions are 
taken to be Gaussian. 

2. {h/e) in a high-density crystal or glass scintillator 

In an EM cascade the electrons are relativistic until their energies fall well 
below the critical energy, so that almost all of the energy is deposited by near- 
minimum ionizing electrons. No appreciable energy exits from the EM cascade 
via photonuclcar interactions. The result is a response very nearly linear in the 
incident electron or photon energy. 

Hadronic interactions deposit energy in a variety of ways. (An inventory is 
given in Table 1 [by Gabriel and Schmidt] in Ref. [9|, and detailed discussions 
can be found in Rcfs. [2i"26| and other recent reviews). A large fraction of the 



hadronic energy (sa 20% for Fe/scintillator and « 40% for U/scintillator sam- 
pling calorimeters) goes to nuclear dissociation and recoil, and is "invisible." 
Neutrinos and most muons escape. Some fraction of the neutrons can be de- 
tected via n-p scattering in hydrogeneous materials such as organic scintillator, 
but much or most of the neutron energy is also lost. Low-energy protons and 
charged fission fragments produce saturated signals in scintillator. (This occurs 



in inorganic l27| as well as organic scintillators [28|.) All of these factors result 
in low visible response to the hadronic component of the cascade relative to 
response to the EM component. 

Detection of recoil protons in neutron scattering in hydrogenous detectors 
increases h [2^. In a sampling calorimeter a disproportionate fraction of the 
EM energy is deposited in the higher-Z absorber; the absorber/active region 
thickness ratio can be "tuned" to decrease e. Both of these effects increase h/e. 
In practical sampling calorimeters (h/e) is typically 0.7, and can be made to 
approach unity with careful design. 

Neither mechanism for increasing h/e is available to a high-density homoge- 
nous calorimeter. 

As we shall see, the resolution is dependent on the "/i/e contrast," the dif- 
ference between (h/e) for the Gherenkov {{h/e\c)) and scintillation( ( h /e\ s) ) 
readouts. Based on experience with quartz-fiber readout calorimeters [30l432l | . 
{h/e\c) — 0.20-0.250 There are few data concerning {h/e\s) in a homogeneous 
calorimeter, but there is no way to hide EM energy in the absorber and there is 



Prom the data shown in Table 3 of Ref. [SQ 
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Figure 1: Energy-independent event locus in the C / E-S / E plane. As indicated by the MC 
events, resolution improves and the mean moves upwrard along the locus as the beam energy 
is increased. 

very little neutron sensitivity. We might expect as much as 30% of the hadron 
energy to be expended on nuclear dissociation and therefore invisible, and 15%- 
20% to be carried by neutrons. These alone would result in {h/e\s) ~ 0.5. There 
are other effects, such as incomplete Chercnkov-scintillator separation and sat- 
urated scintillator response to highly ionizing particles, so (h/e\s) = 0.35-0.5 
might be expected. This is corroborated by a comment in Ref. [16|: "The e/h 
value of ECAL [PbW04] as a scintillation device is much larger than for the 
Cu/plastic sampling structure in DREAM: 2.4 vs. 1.3." (h/e = 0.43 vs 0.7.) 




3. Dual-readout hadronic calorimetry 



In 1997 Wigmans discussed the advantages of adding a quartz-fiber readout 



to the scintillator readout of a sampling calorimeter [33[. With a Cherenkov 
readout (C) fairly blind to hadronic activity and a scintillator readout (S) with 
optimized hadronic response, fem could be determined (or eliminated) and a 
corrected energy found. The DREAM collaboration elegantly implemented this 
proposal with a test-beam calorimeter having quartz and plastic scintillator 
fibers in copper tubes 31,32, [13] . 



In the dual-readout case Eq. [T] is replaced by [Sll, 1351437 1 



S = ElU 
C = E[U 



+ {l-fe,n){h/e\s)] 
+ {1- fem){h/ec)] ■ 



(2) 
(3) 
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In parametric form, Eqns. ([5]) and ([3]) describe a straight line-segment event 
locus in the C-S (or C/E-S/E) plane, as illustrated in Fig. [T] If the cascade 
is "all electomagnetic" (/em — 1), then S/E = C/E = 1. If the cascade is "all 
hadronic" {fem = 0), then S/E = h/e\s and C/E ^ h/e\c. (The MC "events" 
shown in the figure are discussed below.) The slope in the C/E-S/E plane is 
independent of energy; with increasing energy the distribution just moves up 
along the locus. 

It is convenient to introduce the less cumbersome notation h/e\x = rjx- 

S = E[U„ +{1- fern)Vs] (4) 
C = S[/em + (l-/e™)??c] (5) 

Equations and (O, linear in 1/E and /em, can be rewritten as 



S -il-Vs)\fl/E\^fr,s 

C -{I-VC)) {fem J \VC 



(6) 



with solutions 137 



^ ^ ^(1 - r,c) - C{1 - r,s) 
Vs - Vc 

^^"^ - S{l-vc)-C{l-vs) • ^ ^ 

There is an important difference between Eqs. [1] [21 and HI and Eq. [T] The 
first three give estimators of the scintillator response, given /em, VSj and the 
incident energy E. In contrast, Eq. [7] provides an estimator of this energy given 
S, C, rjs and rye- Similarly, Eq. [8] provides an estimator of fem- 

The sensitivity of the energy estimator to the h/e contrast is particularly 
manifest in Eq.[7l If (rys) — (rye) is small compared to the statistical fluctuations 
of r]s and rjc, then the scatter in the energy estimators will be large. 

In Eqs. [7] and [U rjc an 775 are the values peculiar to that event. These are 
unknown — and unknowable, until some way of tagging the hadronic composition 
becomes available!! But in an experimental situation, an estimator of the energy 
must be established for each event. There is little choice but to replace these 
quantities by their means. 

In this case, it is convenient to write the energy estimator (Eq. [7]) more 



compactly as [35l-l37| 



where I have made use of the slope of the event locus (the ratio of ranges of S 
and C) shown in Fig. [T] 



^Neutron detection has been proposed and is being explored for this purpose |38l. l39 
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Figure 2: The mean, standard deviation, and fractional standard deviation of fem in lead, 
as simulated by FLUKA90. This is basically Fig. 3(b) in Ref. [s^. An error in the fit to 
the fractional standard deviation has been corrected, a slightly different MC data set has 
been used, and Wigmans' fit to the fractional standard deviation in a copper/quartz-fiber 
calorimeter has been added [SSl . 



4. Electromagnetic fraction 

In a dual-readout calorimeter estimators of E and /em can be determined on 
an event- by- event basis. The EM fraction fem is no longer a stochastic variable. 

The fem p.d.f. has an energy-dependent mean (near 0.5 at 100-150 GeV, 
approaching unity as £^ — )- oo). Its standard deviation is about 11%, depending 
weakly on energy and calorimeter composition, and it is skewed to the large- /g. 



side 2J, [35| . The mean, fractional standard deviation, and standard deviation 
of fem as simulated by FLUKAQCH for cascades in a large lead cylinder are 
shown in Fig. [2l For this case I obtain the empirical fits 

/ ^ X 0.87-1 
\ -0.076 



(The functional form of Eg. [TT] satisfies the requirement (fem) — > 1 as — ?> oo, 
and a must slowly approach as — > oo.) Although I use these expressions 
for the present calculations, the numerical constants will be somewhat different 
for different calorimeter structures. The fractional standard deviation found for 



*SincG FLUKA90 many improvements have been made, especially in the nuclear physics 
modeling. The high-energy cascade modeling is nearly the same. Only the iz'^ energy fraction 
in the cascade is of interest here. 
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a copper/quartz-fiber calorimeter (Fig. 4.46 in Ref. is shown in Fig. [2]for 
comparison. 

Values of fem must be chosen from distributions with the means and r.m.s. 
widths given by Eqs.[TT]and[T21 and with skewness that agrees with the detailed 
hadronic cascade simulations. This dimensionless skewness 71 (= Hs/cr^, where 
fj,3 is the third moment about the mean) is not well-determined from our simu- 
lations, but it is about 0.6, which I assume here. Any smooth positive function 
bounded between and 1 that can be adjusted to have these properties would 
be satisfactory. I have found it convenient to use the Beta function 

B{x;a,/3) (xx"^\l-xf'\ (13) 

with a and /3 adjusted to obtain the desired a and 71. It is then displaced 
so that its mean is {fem)- The displaced function does not quite go to zero 
at X = 1, but it is sufficiently close to zero that the problem can be ignored. 
The procedure is illustrated in Fig. |3Ua) for E = 200 GeV pions in lead. The 
distribution from the FLUKA simulation is superimposed. 

Values of fem are chosen from this function by choosing uniformly distributed 
random points (/em, y) and retaining the values of fem where y is not above the 
p.d.f. at that value of fem- A histogram of one such array of values is overlaid 
on the model p.d.f. in Fig. [HJ^b). 

Experimental distributions of fem obtained with the DREAM detector are 
shown in Refs. [sij, [32j, and [s^. Because of resolution effects they are consid- 
erably broader and less skewed than the FLUKA-generated distributions. The 
typical fem estimator distribution shown in Fig. ^h) illustrates the resolution 
loss. 



5. Resolution contributions 

In a dual-readout calorimeter fem has been elevated from a stochastic quan- 
tity to a measured quantity. In the absence of other resolution contributions, 
C and S are completely correlated, and the reconstructed energy distribution 
given by Eq. ([7]) or Eq. (jH) is a delta function. 

The signal distribution is broadened by photoelectron (p.e.) statistics, un- 
certain shower leakage corrections, uncorrected signal collection irregularities, 
electronic noise, sampling fluctuations (in the case of a sampling calorimeter), 
intrinsic fluctuations in the visible fraction of hadronic energy deposition, and 
other effects. As a matter of convenience in this discussion, I ignore the contribu- 
tions from the readout-associated factors or assume that their contributions are 
lumped into ctP '^'. Although published resolution measurements often involve 
constants and other deviations from 1/\/E scaling, for our present purposes I 
use 

ffE crP°- cr'"'^ O-'^^P , „ 

— = ^ © ^ © , (14) 

E ^/e Ve Ve 

where the ct's on the right side are fractional resolutions at 1 GeV if E is 
in GeV. The hadronic intrinsic and sampling contributions (ct™*' and a^^™^) 
are discussed below. 
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Figure 3: Approximation of the fem p.d.f. by a displaced Beta p.d.f. for 200 GeV pions incident 
on a lead "calorimeter." (a) The Beta p.d.f. (dashed curve) with a from Eg. I12l and 71 = 0.6, 
and the distribution displaced to the mean given by Eg. 1111 (solid curve). The superimposed 
histogram is from a FLUKA90 simulation, (b) A histogram of 20 000 MC values chosen from 
this p.d.f., compared with the displaced Beta function. The dashed histogram shows the 
corresponding fem estimator distribution obtained via Eq. |8] for the case (ric) = 0.25 and 
= 0.50. It is broader and less skewed than the input distribution. 
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The nuclear dissociation energy deposit distribution is distinctly non-Gaussian; 
for example, see Fig. 4 in Ref. [33| . But visible energy deposit via ionization by 
charged particles is more important in materials of interest, and its distribution 
is not completely correlated with that of nuclear dissociation and other miss- 
ing energy. The resulting visible energy deposit distribution in compensated 



sampling calorimeters seems to be near-Gaussian, and Drews et al. |12| as- 
sume Gaussian distributions in separating intrinsic and sampling contributions. 
Gaussian distributions are assumed in Eg. 1141 and elsewhere in this analysis. 

1. Readout statistics. Photomultipliers or avalanche photodiodes will likely 
be used to detect the scintillation and Cherenkov light. In this case Gaussian 
(Poisson) fluctuations with variance equal to the number of detected photons, 
and hence proportional to the energy, are added to the signals. The fractional 
uncertainty thus scales as 1/Ve. 

Electromagnetic calorimeters provide guidance about obtainable resolution. 
Some of these are tabulated in Sec. 28.9.1 of Ref. [s^. Although fractional reso- 
lution as good as 2%/^/E has been obtained in scintillating crystal calorimeters, 
the best reported Cherenkov response in lead-glass EM calorimeters was 5%, for 
the OPAL endcap [i^, corresponding to 400 p.e.'s/GeV. To maximize the num- 
ber of photoelectrons collected, w 45% of the ends of the lead-glass blocks were 
covered by the PMT's. Such collection efhciency will probably not be possible 
with the proposed crystal calorimeter. 

At the other extremMhe DREAM detector obtained 35%/ VE (8 p.e./GeV) 



with quartz fibers [41|-|43[. Improving Cherenkov light yield has been a major 
goal of the collaboration's recent work with a variety of crystals Il 5|-|22l| . Yields 
as high as 55 p.e./GeV have been reported for Mo-doped PbW04 [21|. Further 
improvements are possible, for example by extending the UV response of the 
photodetectors and UV transmission of the radiators. 

In this model study it is realistic to assume a middle ground, (J^°' = 
o'o.c/^' where (7^;^ ~ 10-20%, or 100-25 p.e.'s/GeV. In most of the MC's 
reported here, 15%, or 44 p.e.'s/GeV was used. 

In most of the dual-readout homogeneous calorimeters under discussion, 
the Cherenkov and scintillator signals are to be separated by color and pulse 
shape. This requires that the scintillation signal is not large compared with the 
Cherenkov signal, say (S) = (C) , where is "a few." We take = 4 in this 
study: c7^-g = 

2. Sampling and intrinsic fluctuations. In a sampling calorimeter most of the 
energy is deposited in the absorber, and there arc large fiuctuations in the 
small fraction of the visible hadronic energy deposited in the active medium. A 
homogeneous calorimeter is not subject to these sampling fluctuations; in fact, 
this is an important reason for choosing itl^j 



^Somc suggested schemes, sueh as alternating lead glass Cherenkov planes with heavy glass 
scintillator planes |44| . are totally active but only quasi-homogeneous. In these cases sampling 
fluctuations have some importance. 
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Table 1: Examples of near-compensating sampling hadron calorimeters. For our present 
purposes some calorimeter structure variation and constant terms in the fitted resolution 
have been ignored. 



Calorimeter 


Passive 


Active 


Resolution 


h/e 


Reference 


(Akesson et al.) 


Cu, U/Cu, U 


Scint (2.5mm) 


36%/VS 


0.90 


[45] 


HELIOS 


U (3 mm) 


Scint (2.5 mm) 


34%//B 


0.984 ±0.006 


[46] 


ZEUS FCAL 


U (3.0/3.2 mm) 


Scint (2,5/3.0 mm) 


35%/^/:B 


1.03 


[47,48] 


WA80 


U (3 mm) 


Scint (3 mm) 


67%/V^ 


0.89 


[49] 


(Drews, et al.) 


Pb (10 mm) 


Scint (2.5 mm) 


U%/^ 


0.90 ± 0.01 


[12, 50] 


(Drews, et al.) 


U (3.2 mm) 


Scint (3.0 mm) 


36%/Ve 


0.99 ±0.01 


[12] 


SPACAL 


Pb (4x scint vol) 


1 mm scint fibers 


30%/Ve 


0.87 


[51, 52] 


PCAL* 


Pbt(10 mm) 


Scint (3 mm) 


32%/Ve 


0.89 


[53] 



* E < 6.8 GcV. 



t Every 6th plate is 16 mm thick Fe. 



Much of the hadronic energy deposit is invisible, going to nuclear disasso- 
ciation, the production of unseen neutrons, etc., with consequent "intrinsic" 
fluctuations in the visible signal even if there is no absorber. In all but a few 
dedicated test-beam experiments, sampling and intrinsic fluctuations arc inex- 



tricable. Drews et al. [12[ studied the problem using compensated sandwich 



calorimeters with separate readouts for odd- and even- numbered layers. Since 
sampling fluctuations from layer to layer are independent, the sum of the odd- 
and even-layer signals, and their difference have the same variance. Signals from 
the intrinsic fluctuations are correlated between layers, so sums and differences 
could be used to separate sampling and intrinsic variances. (ctP °- sa 
did not significantly broaden the responses.) In the case of lead plates, the 
sampling contribution was (41.2±0.9)%/-\/E and the intrinsic contribution was 
(13.4±4.7)%/v'£'. In the case of uranium plates, the sampling contribution was 
(31.1 ± 0.9)%//B and the intrinsic contribution was (20.4 ± 2.4)%//E. Lead 
is probably closest to compositions likely to be considered for the homogeneous 
calorimeter, so this result is relevant to the present discussion. 

Since fem fluctuations do not contribute in a compensating calorimeter (by 
definition of "compensating"), its resolution is the result of p.e., sampling, and 
intrinsic variations!^ Examples of near-compensating calorimeter resolution are 
shown in Table [T] We may take ge = 35%/V^ as a representative "best case." 
Drew et al.'s Pb-plate calorimeter achieved 44%/-\/E (italicized line in Table [1]). 



^This is not quite true. Suppose, for example, that {h/e) = 1, and we select a subset of 
events for which h/e fluctuates to 10% below the mean, to 0.9. For these events S is sensitive 
to fluctuations of fem- The distribution oi h/e about unity thus introduces sensitivity to fem- 
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Figure 4: Input p.d.f.'s and MC realizations of the h/e distributions at 75 GeV for {r]c) = 0.25 
and (ris) = 0.40. cr^if' and o-jSt"-' are from Table[2] 



Scaling to this resolution, erg'"'- « (35/44) x 13.4%//E ss We adopt 

this as the fiducial scintillator intrinsic resolution for this study. If (rjs) ^ 1, 
then crjf*''- =0.11 {r]s) /VE. 

The standard deviation a^^^' of rjc is more problematic. Light produced by 
relativistic pions, with smaller contributions from other hadrons and electrons 
produced via nuclear 7-ray interactions, contribute a hadronic component to 
the Cherenkov response. The QFCAL ^ and DREAM [sSl collaborations 
found (77c) = 0.20-0.25. Since most of the scintiUation signal is produced by 
nonrelativistic ionizing particles, 77s and rjc are nearly independent. 

Typically about 35 (relativistic) tt^'s are produced by a 100 GeV pion show- 
ering in Pb|3 The fractional standard deviation is l/\/35 = 0.17 at 100 GeV, 
or 1.7/Ve. This broad distribution is scaled down by (ijc) ~ 0.20-0.25: 
ab"*-«1.7(77c)/Vl. 

These results are summarized in Table [21 Typical distributions of ?/s and r\c 
are shown in Fig. |4l Since the energy distributions reported for compensating 
calorimeter, e.g., as given in Table [1] are all consistent with Gaussian distribu- 
tion, it is evidently valid to consider the distributions discussed in this section as 
Gaussian as well. This results in near-Gaussian energy-estimator distributions 
as reconstructed from S and C. 



Section 2.3.1.3 of Ref. \2A 
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Tabic 2: Summary of resolution contributions. Reference values in the last column are used 
in Sec. |6] 



Std dcv Value at E Reference value Optimistic value 



p.c. 






p.c. 


15% 


p.c. 


5% 


p.c. 






C = 2 




e = io 






<c- ivc) /^ 


/E 


<c = 


170% 




:85% 




<s- ivs) /x 


fE 


<§■ = 


11% 


5.5% 





6. Simulation of the energy estimator distribution 

Equation [5] and the equivalent Eq. [7] describe the estimator of incident pion 
energy as obtained from the observed scintillation and Cherenkov signals. We 
incorporate the resolution contributions discussed in Sec. [5] and summarized in 
Table[2]to obtain distributions of S, C, and the energy estimator E via a simple, 
transparent Monte Carlo calculation as follows: 

1. Choose the incident energy E and the detection efficiency ratios (775) and 
(ijc)- fixing the latter in the range 0.20-0.25. 

2. Choose the resolution parameters, using the reference values given in Ta- 
bled 

3. Generate an array of N values of /em chosen from the displaced Beta 
distribution for energy E. 

4. Generate N values of 775 from a normal distribution with mean (775) and 
fractional standard deviation cr;^''^'. Similarly, generate values of 77c 
from a normal distribution with mean (rjc) and standard deviation cr™'"''. 

5. Construct the corresponding S and C arrays via Eqs. [4] and[5] 

6. Replace each S and C as calculated in step 4 with values chosen from nor- 
mal distributions with means S and C and standard deviations ctq ^ 
and a^'g^/S, respectively. The resulting S and C include p.e. statis- 
tics, and are used for subsequent "data analysis." (Since aQ ^/VC and 

a^'°g/\fS are fractional standard deviations, a^''(j\fC and Cg'sV^ have 
the same units as S and C (GeV).) 

7. Find the energy estimator array via Eq. [7]or Eq. HI 

The results of four 10 000 event simulations are shown in Fig. [SJ at 75 and 
200 GeV using realistic (0.4) and optimistic (0.6) values of 775. In each case the 
mean value of the estimator of E (E^st) scaled by the beam energy is 1.00. The 
distributions agree well with the Gaussians with the same mean and standard 
deviation drawn over the E histograms. In nearly all cases |7i| < 0.05. 

The fractional standard deviation of the energy estimator scales as 1/ VE. 
The coefficient (resolution at 1 GeV) as a function of h/e\s (= 77s) is shown in 
Fig. |6]for two values of h/e\c (= ijc) that bracket the range found for quartz- 
fiber readout calorimeters I30l - l32t . 
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Figure 5: Monte Carlo distributions of C, S, and the estimator of E. (a) and (c) are for beam 
energies of 75 GeV; for (b) and (d) the beam energy is 200 GeV. For (a) and (b) 775 = 0.40, 
while for (c) and (d) r]g = 0.60. In all cases Tjc = 0.25. Gaussians with the "measured" 
cr^ and mean relative to beam energy are shown as dotted black lines. Details are given in 
Tabic [3] 
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Table 3: Parameters for the examples shown in Fig. [5] eaeh based on 20 000 simulated events. 
In all eases rjc = 0.25. 





(a) 


(b) 


(c) 


(d) 


Eb 


75 GcV 


200 GcV 


75 GeV 


200 GcV 


Vs 


0.40 


0.40 


0.60 


0.60 


(S/Eb) 


0.67 


0.71 


0.78 


0.81 




11.6% 


10.1% 


6.7% 


5.9% 


71 s 


0.57 


0.62 


0.58 


0.59 


(C/Eb) 


0.59 


0.64 


0.59 


0.64 


oc 


17.2% 


14.2% 


17.2% 


14.2% 


7ic 


0.42 


0.55 


0.45 


0.53 


{Ecst/Eb) 


0.999 


1.000 


1.000 


1.000 




12.8% 


7.24% 


3.96% 


2.27% 


as @ 1 GeV 


111.% 


102.% 


34.3% 


32.1% 


JlE 


0.008 


0.002 


-0.015 


-0.003 



7. Discussion and conclusions 

The standard deviations described in this section are highly uncertain, and 
can be used only as guides. The effect of large changes in the input variances on 
the resolution curves shown in Fig. [6] have been studied. Examples of optimistic 
excursions from best estimates are given in Table[2] At (775) ~ 0.45 and (rjc) = 
0.20; where ctb = 52.5%/ Ve, the following improvements are found: 

1. Increase the Chercnkov p.e. yield by a factor of nine, from cTq c; ~ 

(44 p.e.'s/GcV) to crP °- = 5% (400 p.c.'s/GeV). The resolution improves 
from 52.5%/ VE to 47.1%/v^, or -5.4% — not a large improvement. 

2. Increase the scintillator p.e. yield by a factor of 25 (^ = 2 — > ^ = 10). The 
resolution improves from 52.5%/a/E' to 49.0%/v^, or -3.5%. At least in 
this model, there is no serious penalty for using a "weak" scintillator. 

3. Decrease the width of the Cherenkov intrinsic hadronic resolution at 1 GeV 
by a factor of two, from 170% to 85%. The former number is expected from 
fluctuations in the number of relativistic pions produced in the cascade, 
and so should be relatively dependable. However, this change produces a 
large improvement: 52.5%o/VE — > 39.1%o/VE, or —13.4%. The curve for 
this case has been added to Fig. [6] 

4. Finally, halve the width of the scintillator intrinsic hadronic resolution 
fromll%/v^to5.5%//E. As expected, there is little effect: 52.5%/V^^. 
52.1%, or —0.4%. A substantial downward excursion of this number would 
be unphysical. 

The overall conclusion remains and is evident from Eq. [T] in a homogeneous 
hadronic calorimeter with dual readouts and a non-hydrogenous scintillator. 
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Figure 6: Dependence of resolution on rjs at two values of rjc = h/e\c- The results for 
r]c = 0.20 are essentially those for rjc = 0.25 displaced to the left by 0.05; it is the contrast 
between rjQ and rjg that determines the resolution. The dashed line shows the effect of halving 
the width of the riQ p.d.f. 

event-to-event fluctuations of rjc and rjs destroy the resolution if their means 
are close. A large contrast between them is unlikely because the mechanisms 
that increase (77s) are not present. It remains true that the mean of the en- 
ergy estimator distribution is the beam energy, and its distribution is nearly 
Gaussian. 

Detailed simulations with a modern, sophisticated Monte Carlo program 
such as GEANT4 would be desirable before commitment to a large R&D effort. 
It would be particularly informative to understand the means and variances of 
77c and 77s. 
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